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Its  abstract  .  1 

A  procedure  for  determining  the  mean  and  covariance  errors  in  an  oh?  filter  operating 
in  cartesian  coordinates  was  found.  Hie  results  obtained  from  this  procedure  were  com¬ 
pared  to  an  filter  operating  in  polar  coordinates. 

Assuming  that  the  Input  measurements  in  polar  coordinates  were  Gaussian  distributed, 
it  was  drown  that  at.  the  output  of  the  coordinate  transformations  the  noise  could  be  ap¬ 
proximated  accurately  by  a  Gaussian  distribution  for  typical  radar  data.  Oo6ed-farm  solu¬ 
tions  under  steady -state  conditions  were  found  for  the  output  covariances  for  the  polar 
coordinate  filter  and  for  the  cartesian  coordinate  filter  when  the  target  is  stationary.  These 
covariances  depended  upon  a,  0,  and  the  measurement  variances.  For  moving  targets,  the 
cartesian  coordinate  filter  yielded  output  covariances  which  were  nonstationary.  Their 
values  depended  upon  a,  $,  measurement  variances,  target  trajectory,  target  speed,  and 
sampling  time.  The  mean  errors  were  discussed.  Under  fading  conditions  both  the  mean 
and  covariance  errors  increased  during  the  fading  time. 
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ABSTRACT 


A  procedure  for  determining  the  m end  covariance  errors  in 
an  a-0  filter  operating  in  cartesian  coordinates  was  found.  The  re¬ 
sults  obtained  from  this  procedure  were  compared  to  an  a-0  filter 
operating  in  polar  coordinates. 

Assuming  that  the  input  measurements  in  polar  coordinates  were 
Gaussian  distributed,  it  was  shown  that  at  the  output  of  the  coordi¬ 
nate  transformations  the  noise  could  be  approximated  accurately  by 
a  Gaussian  distribution  for  typical  radar  data.  Closed-form  solutions 
under  steady-state  conditions  ware  found  for  the  output  covariances 
for  the  polar  coordinate  filter  and  for  the  cartesian  coordinate  filter 
when  the  target  is  stationary.  These  covariances  depended  upon  a, 
i J ,  and  the  measurement  variances.  For  moving  targets,  the  cartesian 
coordinate  filter  yielded  output  covariances  which  were  nonstation¬ 
ary.  Their  values  depended  upon  a,  0,  measurement  variances,  target 
trajectory,  target  speed,  and  sampling  time.  The  mean  errors  were 
discussed.  Under  fading  conditions  both  the  mean  and  covariance 
errors  increased  during  the  fading  time. 
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DESCRIPTION  OF  AN  ct-0  FILTER 
IN  CARTESIAN  COORDINATES 


1.0  INTRODUCTION 

In  the  last  several  years,  there  has  been  a  considerable  amount  of  interest  in  auto¬ 
matic  detection  and  tracking  for  search  radar  systems.  Several  systems  exist  with  varying 
degrees  of  automation  such  as  MTDS,  NTDS,  and  the  S PS-33.  Others  are  being  proposed 
such  as  the  Gillfillan  and  APL  systems  for  the  SPS-48,  the  JPTDS  program  for  the  SPS-49, 
and  the  AEGIS  system.  Even  with  these  efforts  there  is  still  a  need  to  improve  system 
performance  under  various  conditions. 

NRL  Report  7434  recently  studied  the  effects  of  maneuvering  targets,  measurement 
noise,  false  targets,  and  fade  conditions  on  the  ability  of  an  a-0  filter  operating  in  polar 
coordinates  to  maintain  a  track  (1).  Even  mere  recently  NRL  Report  7505  discussed  the 
ability  of  this  polar  coordinate  filter  to  hand  off  its  track  from  the  search  radar  to  the 
track  radar  (2).  In  both  of  these  reports  a  considerable  amount  of  difficulty  was  en¬ 
countered  in  either  maintaining  or  handing  off  a  track  at  close  ranges  when  a  polar  co¬ 
ordinate  filter  was  used.  This  was  due  to  large  range  and  azimuth  accelerations  at  the 
dose  ranges.  In  the  cartesian  coordinate  system  these  large  accelerations  are  not  en¬ 
countered.  However,  the  nonlinear  tra.<&£ormations  encountered  between  the  two  coor¬ 
dinates  change  the  noise  processes,  lit  is  the  purpose  of  this  report  to  describe  analytically 
the  a-0  filter  operating  in  cartesian  coordinates  and  compare  these  results  with  the  results 
of  the  polar  coordinate  filter. 

Section  2.0  describes  the  probability  densities  under  the  coordinate  system  transfor¬ 
mations.  Section  3.0  describes  the  general  characteristics  of  the  filter  and  the  mean 
errors  between  the  predicted  and  true  target’s  positions.  Section  4.0  describes  the  covari¬ 
ances  at  the  output  of  the  filter  system.  Section  5.0  studies  *he  mean  and  covariance 
errors  under  fading  conditions  and  presents  the  results  of  a  simulation  calculating  the 
probability  of  placing  tne  beam  of  the  tracking  radar  on  a  target  using  the  track  set  up 
by  the  <x-0  filter.  Conclusions  are  given  in  Section  6.0. 


2.0  THE  NOISE  PROCESS 

In  the  study  of  any  filter  it  is  essential  to  know  the  characteristics  of  the  desired 
signals  and  the  noise  which  excite  the  filter.  The  mean  motion  of  the  targets  is  studied 
in  Section  3.0.  The  description  of  the  noise  processes  proceeds  as  follows: 

The  block  diagram  of  the  filtering  system  is  shown  in  Fig.  2.1.  The  polar  coordinate 
radar  measurements  are  Rm  in  range  and  6,~  in  azimuth,  where  Rm  and  6m  are  assumed 
to  be  uncorrelated,  Gaussian,  amplitude-distributed  random  variables  with  means  6m,  Rm 
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fig.  2.1— FUter  system 


and  variances  aBm>  °lm*  In  addition,  the  measurements  are  assumed  to  be  independent 
from  scan  to  scan  of  the  search  radar.  ThL  section  is  concerned  with  determining  ap¬ 
proximate  probability  densities  p(Xm,  Ym),  p(Xp,  Yp),  and  p(Rp,  op). 


1 

?' 


2.1  Polar  to  Cartesian  Coordinate  Transformation 

The  probability  density  of  the  polar  coordinate  radar  measurement  is 


P(Rmt  @m)  ~ 


2*  °Rm  °en 


exp 


f  1 

(Rm  ~  Rmf 

— x 

.  (®m  —  &m) 

p 

°R 

_  Km 

°m  j 

(2.1) 


Contours  describing  constant  values  of  the  probability  density  function  are  plotted  for 
two  different  cases  in  Figs.  2.2  and  2.3.  Observing  the  central  (111*8)  regions  of  these 
densities,  one  finds  that  this  region  appears  to  be  a  correlated  Gaussian  process  in  the  . 
(Xm ,  Ym }  coordinates.  This  observation  is  next  investigated. 

A  cartesian  coordinate  system  (p,  q)  is  defined  as  shown  in  Figs.  2.2  and  2.3.  For 
an  arbitrary  point  (Rm,  Qm)  in  polar  coordinates,  the  values  of  p  and  q  in  the  p-q  rectan¬ 
gular  coordinate  system  are  found  to  be 

P  ~  -®m[2  ~  cos  {d/ft  —  Off.,  j  —  Rm  (2.2) 


q  —  Rm  sin  (hn)  (2.3) 

with  the  aid  of  Fig.  2.4,  For  rases  when  ($*  -  )  is  less  than  about  5°,  Eqs.  (2.2)  and 

(2.-3)  can  be  approximated  by 


P  ~  Rm  Rm 

(2.4) 

«?  =  Rmtfm  ~  §n) 

(2.5) 

with  very  little  error.  For  example,  if  =  0.5  degrees,  one  would  be  at  10  a$m  or  in 
the  far  tail  region  before  tire  approximation  begins  to  be  significantly  in  error.  Further¬ 
more,  if  does  not  significantly  deviate  from  Rm ,  one  can  further  approximate  Eq. 
(2.5)  as 


Q  *"  Rm'^m  &m)  - 


(2.6) 
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q  =  Rm  s M9m-Sn) 


Fig.  2.4— Geometry  required  to  compute  p  and  q 
in  terms  of  polai  coordinates 

For  example,  if  ORm/Rm  -  0.01  as  would  be  the  case  for  Rm  =  4.18  nani.  and  ORm  of 
250  ft,  one  would  be  in  error  by  1%  at  one  standard  deviation  and  5 %  at  five  standard 
deviations.  At  longer  ranges  the  error  is  much  less. 

Equations  (2.4)  and  (2.8)  are  linear  transformations  and  therefore  p  and  q  are 
Gaussian  distributed,  at  least  to  the  extent  in  which  the  approximations  are  valid.  By 
rotating  die  (p,  q)  coordinates  and  shifting  the  mean,  one  obtains  the  (X-.  Y-t) 
coordinates 

Xm  ~  Xm  -  p  cos  0m  —  q  sin  6m  (2.7, 

Ym  -  Ym  ~  P  sin  0m  +  q  cos  6m  .  (2.8 

Again,  these  are  linear  transformations  and  therefore  tfa  i  variables  {Xm,  Ym)  are  Gaussian 
distributed: 
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P&ms  Ym ) 


2,r  °xm 


where 


r  r  _o 
(Xm-Xmf 


X  exp  - 


-  2a 


xmym 


(Xm  ~  Xm)(Ym  +  (Kn-^n)2 


0lm 


1  -  p. 


%ym 


(2.9) 

< 

=  °Rm  cos20«  sin20m 

(2.10) 

°L 

'  °Rm  sin2®m  +  cos2^ 

(2.11) 

Pxmym 

[°lm  “  (^V)2]  S&n2&m 

2  ®xm  °ym 

(2.12) 

xm 

=  Rm  cos  0-„ 

(2.13) 

Ym 

=  Rm  sin  0m. 

(2.14) 

An  independent  procedure  for  obtaining  Eqs.  (2.10)  through  (2.14)  is  also  shown  in  the 
appendix. 


The  two  density  functions  used  in  plotting  Fi^.  2.2  and  2.3  are  approximated  by 
p(Xm,  Ym)  shown  in  Eq.  (2.9),  and  .he  results  are  shown  in  Figs.  2.5  and  2.6.  In 
comparing  Fig.  2.2  with  2.5  and  2.3  with  2.6,  one  finds  that  the  central  and  near  tail 
regions  of  the  two  densities  p{Xm,  Ym)  and  p(Rm,  0m)  are  essentially  the  same.  The  ap¬ 
proximation  begins  to  break  down  in  the  far  tail  regions.  However,  it  is  the  central  and 
near  tail  regions  which  control  the  system.  Any  point  in  the  far  tail  usually  results  in 
saturation  or  in  this  case  no  target  being  accepted. 

2.2  Effect  of  Linear  Filter 

Since  p(X, Ym)  is  essentially  Gaussian,  the  output  of  a  linear  filter  is  also  Gaussian 
distributed;  i.e.. 
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P(XP,  Yp)  = - =—=== 

2lraxpOypy/^~Pxpyp 


h  (2-15) 


This  is  because  the  output  of  any  filter  can  be  written  as  a  linear  ccmbination  of  the 
inputs,  and,  since  the  sum  of  Gaussian  random  variables  is  Gaussian  distributed,  the  out¬ 
put  of  the  filter  is  Gaussian.  The  means  and  covariances  will  be  investigated  in  detail 
later. 

2.3  Cartesian  tc  Polar  Coordinate  Transformation 

Contours  describing  constant  values  of  the  probability  density  function  given  in  Eq. 
(2.15)  are  plotted  for  two  different  cases  in  Figs.  2.7  and  2.8.  Observing  these  figures,  it 
appears  that  the  central  and  near  tail  regions  of  these  densities  are  Gaussian  distributed  in 
polar  coordinates.  This  conjecture  is  investigated.  The  (p,  q)  axis  system  is  defined  in 
Figs.  2.7  and  2.8.  The  {Xp  -  Xp),  (Yp  -  ^p)  axis  is  rotated  so  as  to  coincide  with  the  p-q 
axis.  Then, 


'VV 


X  exp 


(• r 


■2 


’xpyp 


(Xp-Xp)(Yp-Yp)  (Yp  -  Yp)2 

/V  /T  O 


*p  yp 


JyP 


i  -  p 


xpyp 


Fig.  2.7— Constant  contours  of  p(Xp,  Yp) 
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Pig.  2.8-Constanc  contours  or p(Xp,  Yp) 


p  =  {Xp  -  Xp)  cos  6p  +  (Yp  -  Yp)  sin  Qp  (2.16) 

q  =  ~(Xp  -  Xp)  sin  dp  +  (Yp  -  Yp)  cos  6p  .  (2.17) 

For  small  angular  and  range  deviations  in  the  coordinates  (0P,  Rp),  the  approximations  de¬ 
fined  in  Eq.  (2.4)  and  (2.5)  are  valid: 


£ 

1 

■s01 

II 

•a 

(2.18) 

Rp(dp  Op)  -■  q  . 

(2.19) 

Combining  Eqs.  (2.16H2.19),  one  obtains 

Rp  -  Xp  cos  Op  +  Yp  sin  dp 

(2.20) 

dp  -■  (~Xp  sin  dp  +  Ip  cos  dp)/Rp  . 

(2.21) 

Since  Eqs.  (2.20)  and  (2.21)  are  linear  transformations  on  Gaussian-dismbuted  random 
variables  Xp  and  Yp,  p(Rp,  0p)  is  Gaussian  distributed  at  least  over  the  region  in  which  the 
approximation  is  valid: 
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P(Rp-.  dp)  - 


2*%  °Op\/i-Pr 


Rpdp 


where 


X  exp 


f 

(Rp  -  Rp)2  _  (RP-RP)(eP-dp) 

$ 

i 

_ i 

-  1 

o*p  ipVp  oRp  O0p 

of 

Op 

2 

L  1 "  ph*P 

r 

✓ 

,  (2.22) 


aRp  =  axF  cos2eP  +  2°xp  %  Pxpyp  COS  Qp  sin  Op  +  oj  sin\ 
°xp  sin2^>  "  2 oXp  Oyp  pXpyp  sin  6p  cos  \  +  a2^  cos20p 


°i= 


PPpOp 


% 


«>*{ °Xp  +  <4)  “*  29P  +  Pxpyp  0Xp  Oyp  COS 


2 9„ 


ORp  oep  Rp 


Rp  =  y/x*  +  Y 2 


dp  =  tan-1  Yp/Xp  . 


(2.23) 

(2.24) 

(2.25) 

(2.26) 
(2J27) 


The  two  density  functions  used  in  plotting  Figs.  2.7  and  2.8  are  approximated  by 
PiRp,  dp)  shown  in  Eq.  (2.22),  and  the  results  are  shown  in  Figs.  2.9  and  2.10.  In  com- 
parmg  Fig.  2.7  with  2.9  and  2.8  with  2.10,  one  finds  that  the  centra!  and  near  tail  regions 
of  the  two  densities  p{Xp,  Yp)  and  p(Rp,  6p)  are  essentially  the  same.  Again  as  in  Section 
2.1,  the  approximation  begins  to  break  down  in  the  far  tail  regions.  However,  these 
regions  are  of  little  interest. 


2.4  Discussion  of  Results 

The  section  showed  that  if  p(Rm,  8m)  was  Gaussian  distributed  with  small  variances, 
p{Rp,  dp)  was  also  Gaussian  distributed  over  the  central  and  near  tail  regions  of  the  dis¬ 
tribution.  The  means  and  variances  after  each  of  the  coordinate  transformations  were 
found.  The  means  and  variances  following  the  linear  filter  will  be  investigated  later. 


3.0  FILTER  DESCRIPTION 

The  a-fi  filter  is  defined  in  this  section  and  a  few  of  its  characteristics  are  shown. 


NRL  REPORT  7548 


1 


3.1  Filter  Definition 

The  filter  in  the  x  coordinate  is  described  by 


-m  (1  -  8) 


~X(k)~ 

Vx(k) 

t_ 

”  X{k  -  1  )~ 

a 

_vx(fe  - 1>_ 

PIT 

[Xp(fe+;,^=[1  jTItn] 


X(k} 

Vx(k) 


[^(fe)]  (3.1) 


,  for  j  =  1,  ...  m  .  (3.2) 


Similarly,  the  description  in  the  y  coordinate  is 


~Y{k)~ 

"(l-«) 

(1  -a)T 

T  Y(k  - 1) 

“  — 

a 

Vy  (ft) 

-PIT 

(1-P)  _ 

[><*-!>_ 

_pit_ 

[Yp(k+jlm)]  =  [l  j'T/mj 


Y(k) 

Vy(k) 


[ym(fe)3  (3.3) 


for  j  =  1,  ...  m  .  (3.4) 


Since  the  equations  are  identical  in  each  of  the  coordinates,  it  is  sufficient  to  show  a  few 
of  the  characteristics  in  the  x  coordinate. 


3.2  Frequency  Response 

Using  z-transform  analysis,  we  define  the  transfer  functions  (3)  as 


az 


Gx  = 


Gv  = 


X(z)  _ _ 

Xm(z)  z2  -  z(2-a-f} )  +  (1-a) 

Vs(z)  _  (PIT)  z(z-  1) 

*m(z)  ~  z2  -  z(2~a~P)  +  (1-a)  ’ 

and  for  j  =  m  we  can  define  the  transfer  function 


*p(z) 


(a  +  P)  z 


Gp  *m(z) 


z2  -  z( 2-a-p)  +  (1-a) 


(3.5) 

(3.6) 


(3.7) 


By  placing  z  =  eJb}T  into  Eqs.  (3.5M3.7),  the  magnitude  and  phase,  defined  as 
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magnitude  iGxi  |Gt.| 
phase  4>x  4>v 


\GP\ 

$p  » 


can  be  found  as  a  function  of  a*  0,  and  coT.  Gx  and  Gv  are  plotted  for  a  given  case  in 
Fig.  3.1.  This  figure  shows  that  X(k)  is  the  result  of  passing  A"m(fe)-through  a  low-pass 
filter,  and  Vx(k)  is  the  result  of  differentiating  Xm  (k).  The  frequency  of  the  input  signal, 
the  sampling  time  T,  and  the  filter  parameters  a  and  0  control  the  filter’s  response. 

It  is  useful  to  place  Eq.  (3.5)  into  the  form  of  a  classical  second-order  system  (3): 

■  \  «  ) 


az 


G~  = 


z2  -  z  2  exp  (~£c cqT)  cos  u ijT  +  exp  (-2 £co0T) 


(3.8) 


Fig.  3.1— Frequency  response  of  o  0  filter 

Equating  coefficients  between  Eqs.  (3.5)  and  (3.8),  results  in 

a  =  1  -  e~2?w0 T 

0  =  1+  e~2?w07‘  -  2  cos  cojT  . 


(3.9) 

(3.10) 


The  inverse  relations  are 
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En  (iMTi) 


2-^1  “  “ 


(3.11) 


(3.12) 


and 

c*  =  <od/VT^  13.13) 

where  £,  wd.  and  a>o>  are  the  classic  damping  coefficients,  damped  natural  frequency,  and 
natural  frequency  of  a  second-order  system. 


3.3  Errors  Under  Sinusoidal  Excitation 

A  target  having  a  circular  motion  is  used  to  represent,  a  turning  target.  The  geometry 
is  shown  in  Fig.  3.2.  The  equations  of  motion  are 

Xm  =  -£n  +  \Xm\  cos  co0f  <3-14> 


and 

Ym  =  Ym  *  l^nl  sin  WQt  .  (3-15) 


Fig.  3.2— Target  geometry 
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The  waveforms  produced  by  the  circular-motion  target  are  passed  through  the  filter  de¬ 
scribed  by  Eqs.  (3.1  H  3.4).  Since  Am,  Ym  and  the  operations  are  well  defined,  it  is  pos¬ 
sible  to  obtain  a  closed-forrn  solution  for  R(k  +  Urn)  and  9(k  +  jjm).  However,  the 
closed-form  solutions  are  lengthy  and  involved.  It  is  easier  to  simply  compute  numerically 
the  results  under  various  conditions.  The  error  between  the  predicted  position  Xp(k  +  j/m) 
and  the  true  target  position  is  computed  as  a  function  of  time  as  shown  in  Fig.  3.3.  The 
envelope  of  the  peak  error  is  sinusoidal,  as  would  be  predicted  from  ‘inear  system  theory. 
In  addition,  the  envelope  of  the  lower  peaks  is  almost  sinusoidal  and  would  be  if  j  =  0. 

The  error  shown  is  valid  only  a\  the  sample  instants. 


I2f  g  =  96  6  ft.'s2 

v  -  !!00  f»/s 
a  -  05 


Fig.  3.3— Error  between  predicted  and  true  position  for  a  circular-motion  target 


The  errors  between  the  predicted  and  true  ranges  and  the  predicted  and  true  azimuths 
are  computed  next.  Various  examples  are  shown  in  Figs.  3.4,  3.5,  and  3.6.  Ir,  these  fig¬ 
ures,  we  find  that  the  amount  of  filtering  affects  the  envelope  of  the  peak  error  and  the 
ripple  errors.  Although  not  shown,  the  target  trajectory7  and  the  sampling  tine  also  affect 
the  error. 

3.4  Comparison  of  Mean  Errors  for  Cartesian  and  Polar  Coordinate  Filters 
The  a-&  tracker  in  polar  coordinates  is  described  by 

R{k) (1  -  a)  (1  -  or)r]  R(k- 1)"  fa  j 

=  +  j  [*«(*)] 

VR(k)  -p/T  (l-jJ)J  Vfl(fe-l)  J/T\ 


(3.16) 
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mg.  3.5— Mean  error  between  predicted  and  true  position  for  a  circular 
target  trajectory  using  an  a-jS  tracker  in  cartesian  coordinates 


6ik) 

'(1-a) 

a-a)f 

r 

rH 

i 

as, 
t  Sa 

ra" 

= 

1 

+  l  j 

i 

« —  _ 

_~PIT 

a-»_ 

w  — 

1  -Mm 

jwti 

Rp{k+jlm )  =  R(k)  <  {jjm)  T  VRik ) 
Gp(k  + jjm)  =  6(k)  +  (jjm)  T  *&(£) . 


(3.17) 

(318) 

(319) 


Is* 

o 


BEN 


target  trajectory  using  an  Q-$  tracker  in  cartesian  coordinates 


For  the  same  circ  iiar  flight  path  as  shown  in  Fig.  3.2,  the  error  between  the  true  and 
predicted  positions  in  range  and  azimuth  is  computed  and  the  results  are  shown  in  Figs. 
3.7  and  3.8.  As  shown  in  these  figures,  the  errors  in  the  polar  coordinate  tracking  system 
at  the  near  ranges  are  larger  than  the  cartesian  coordinate  ones.  This  is  due  to  the  large 


S 

K 

X 


I 

V- 

% 

< 


Fig.  3.7— Mean  error  between  predicted  and  true  positions  for  z  circular 
target  using  an  cc*£  tracker  in  polar  coordinates 
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v  =  1100  f{/s 


Fig.  3.8— Mean  error  between  predicted  and  true  postions  for  a  circular 
target  trajectory  using  an  d-fl  tracker  in  polar  coordinates 


accelerations  set  up  by  the  trajectory  in  polar  coordinates.  At  the  Tar  ranges  we  find  that 
the  errors  in  ether  tracking  system  are  nearly  the  same.  This  point  can  also  be  illustrated 
by  computing  the  errors  using  both  systems  for  a  constant-velocity,  straight-line  flight 
path.  For  the  cartesian  system  in  steady  state  the  error  is  zero.  Figure  3.9  shows  a 
typical  errcr  sequence  for  tracking  in  the  polar  coordinate  system.  Again,  at  the  far 
ranges  th_  mean  tracking  errors  are  essentially  the  same  for  either  system.  At  short 
ranges  the  mean  error  in  the  polar  system  can  become  quite  large  due  to  the  target  mo¬ 
tion  and  orientation  with  respect  to  he  radar. 


3.5  Discussion  of  Results 


This  section  investigated  the  mean  response  of  the  filter.  It  was  briefly  described 
how  et,  p\  T,  target  trajectory,  and  range  affected  the  error  between  the  predicted  and 
true  positions  of  the  target.  It  was  found  that  a-fi  trackers  operating  either  in  polar  or 
cartesian  coordinates,  were  essentially  equivalent  at  the  far  ranges  in  the  mean  errors. 

The  mean  errors  in  the  polar  coordinate  tracking  system  were  much  larger  than  in  the 
cartesian  coordinate  system  at  short  ranges. 

4.0  RESPONSE  TO  NOISE 

This  section  is  concerned  with  computing  the  covariances  of  the  noise  at  the  output 
of  the  filter  after  the  measurement  noise  described  in  Section  2.0  has  passed  through  the 
filter  described  in  Section  3.C. 
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i 

I 


-Cl 


-15*- 

Fig.  3.9 — Error  between  predicted  and  true  target  positions  at  sampling 
instants  j^/m  for  a  straight-iine  target  trajectory  (tracking  in  polar 
coordinates) 


4.1  Covariance  Equations 

Equations  (3.1)-{3.4)  can  be  written  in  the  form  (4) 


where 


W(fe)  =  A  W(ft-l)  +  r  V(fe), 


W(fe)  = 


A  = 


~X(kf 

a 

°1 

vx(k) 

~Xm{k) 

PIT  0 

V(fe)  = 

r  = 

Y(k) 

0 

a 

Vy(k)_ 

— 

P!T_ 

‘(1  -  a) 

(l-a)T 

0 

0 

— 

'PIT 

(1-0) 

0 

0 

0 

0 

(1  -  a)  (1  -  a)T 

0 

0 

-pit  a- 

P)J 

(4.1) 


The  covariance  equations  can  be  written  as 
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p(fe)  =  a  ?(k  -l)  at  +  r  Q(fe)rT , 


where 


PCft)  = 


pxx(k) 

Wfe> 

Pxy{k) 

%(*) 

-v<*> 

puxvx(k' 

*W*> 

*,*(*> 

PyVx{k) 

Pyy(fe) 

Pyvy(fy 

V(fc) 

PUyVX(k) 

V(fe) 

Pvyvy(k) 

pxvx  =  cov(X,  Vx) ,  etc. 


Q(&)  = 


Qxx(k)  Qxy(k) 

QyX{k)  Qyyik) 


Q^ik)  =  cov(Xm,  Xjjj)  =  o|m 

Q^k)  =  cov(Xm,  Ym)  =  pXmym  aXm  oym 

Qyy(k)  =  co v(Ym,  Ym)  =  a2m  . 

The  cov.aiance  equations  for  the  filter  are 


Pxxik) 

(1-ce)2 

2(1-  a)2T 

(1  -  a)2T2 

pxvx&) 

- 

-0(1  -  «)/ T 

(l-a)(l-20) 

(1  -  a)(l  -  0)T 

JWW 

{PIT)2 

-2P(1  -0)/T 

(1-0)2 

P»(*-D 

a2 

PXVX^  ~  1) 

+ 

00/ r 

_PvXvx(k  ~ 

(0/T)2 

[<?«<*>] 


(4.2) 


(4.3) 
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Pw<*) 


(1  -  a)! 


2(1  -  <x)2T 


-  iv\2t2 


(1-a) 


Pyv(k)  =  -0(1 -a )/T  (l-a)(l-20)  (l-a)(l  -0)T 


-2p(l-0)/T 


(1  ~  (3)2 


Pyy(k-1)  0? 

X  Pyu?,(«-1)  +  Ct&IT  [Qyy(fe)] 

P0^(k~  1)  tf/D2 


^(fe) 


(1-a)2 


*W*> 


(1  -  a)2T  (1  -  a)2T  •  (1  -  a)2T2 


0(l-a)/T  (1-a)  (1-0)  -0(l-u)  (l-a)(l-0)T 


Pv,y(k)  -0(l-a)/T  -P(l-a)  (l-P)d  -«)  (l-a)(l-0)T 


(-0/T)(l-0)  H/r)(l-U) 


(1-0)5 


P„x(fc-l)l  f  a2 


V(fe_1) 

x  + 


pWfc"D  tf/T) 


[«**<*)] 


Forming  the  covariances  of  Eqs.  (3.2)  and  (3.4)  yields 

a2  =  P„(s)  +  2  Ulm)TPXVrlk)  +  U/m)2T2PVrV(k) 


a2  =  P„(fc)  +  2(jlm)TPyt,(k)  +  UI™)2T2Pv„(k) 


yp  yy v  '  wr*—/  *  ‘V'y  '  ' 

=  Vfc)  +  <^>2TV(fe)  +  Vk)3 +  Ulm)2?2PVxVy(k).  (4.8) 

For  stationary  noise  inputs  dosed-form  solutions  can  be  found  for  Eqs.  (4.3H4.5)  by 
pinring  P(fe  +1)  =  P(.%)  and  solving  the  resulting  algebraic  equations.  However,  it  is  easier 
to  obtain  solutions  by  recursively  solving  Eqs.  (4.3)-(4.5)  until  a  steady-state  solution  is 
obtained.  It  is  necessary  to  eliminate  T  as  a  parameter  by  substituting  Ay  =  VyT  and 
Ax  =  VXT  into  the  original  filter  equations.  The  resulting  covariance  equations  are  of  the 
same  form  as  Eqs.  (4.3H4.8)  with  T  =  1.  (The  covariances  are  independent  of  sampling 
time.)  Solving  Eqs.  (4.3)*(4.8)  yields  the  following  result: 


NRL  REPORT  7548 


21 


F  =  °xpIQxx(k)  =  fykVQyyih)  =  %*p%%/G*y(*>  .  (4.9) 

For  all  admissible  a.  /J,  and  j;  for  T  =  1;  and  for  F  -  f  (a,  /3,  ;'). 

A  simple  procedure  for  determining  the  solution  for  stationary  targets  is  next 
described. 


4.2  A  Simple  Stationary  Solution 

For  the  system  shown  in  Fig.  2.1  the  measurement  means  and  covariances  are  trans¬ 
formed  into  ff, rm,  Oym,  Pxmym,  Xm,  Ym  by  Eqs.  (2.10)-(2.14).  One  then  computes  the 
covariances  at  the  output  of  the  filter  by  Eqs.  (4.3)-(4.8).  The  results  a|  ,  Ogp,  pRpdp, 
and  Rp  are  then  computed  by  Eqs.  (2.23)-(2.27).  If  the  target  is  stationary  the  input 
process  to  the  filter  is  stationary  (Eqs.  (2.10)-(2.14)),  and  therefore  the  output  of  the 
filter  has  a  stationary  solution  given  by  Eq.  (4.9).  In  addition,  6m  =  8p  and  =  Rp. 

The  output  of  the  filter  can  then  be  written  as 

°Xp  =  F  °Rm  cos2  8p  +  F(RpaeJ  sin20p  (4.10) 


ayp  =  F  °L  sin2  °p  +  F(Rpaem)2  cos20P 

Pxpyp%°yp  =  05  F  °xm°ym  (°Im  ~(Fp  °oJ2)sin  2®P 
Substituting  Eqs.  (4.10)-(4.12)  into  Eqs.  (2.23)-(2.27)  results  in 


(4.11) 

(4.12) 


°R 

HP 

ii 

3 

(4.13) 

Op 

=  Fag 

“m 

(4.14) 

PRpbp 

=  0 

(4.15) 

RP 

= 

(4.16) 

=  . 

14.17) 

and 


But  this  is  the  same  solution  as  would  have  been  obtained  if  the  system  shown  in  Fig.  4.1 
had  been  used.  The  equations  describing  the  system  are 


~R(k)~ 

yR(k)_ 

-0/r  (l-£) 


~R(k- 1)~ 

a 

+ 

J4(fe-1)_ 

jn 

[*«<*)] 


(4.18) 
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LINEAR 

FILTER 


% 


Fig.  4.1— Equivalent  system  tc  Fig.  2.1  under 
linear  approximations  used  in  Section  2.0 


and 


6(k)  | 

Vfl(fe)j 


<l-«) 

~PIT 


(1~P) 


~9(k-l) 

a 

+ 

Vg(k- 1) 

— 

i 

PIT 

Rp{k  +  j/m )  =  R(k)  +  ( j/m)TVR(k ) 
0p(k  +  j/m)  =  6(k)  +  (jlm)TVe(k). 


The  covariance  equations  are  formed  in  the  same  manner  as  before: 


PRR(k) 

(l-0£)2 

2(1  -  a)2T 

(1  -a)2r2 

PrVrW 

= 

-0(1  -  a)/T 

(1  -  cr)(l  -  2p) 

(1  -  a)(l  -  P)T 

pvRvR(k ) 

(PIT)2 

—2/5(1  -  P)/T 

(1  -P)2 

Prr&-  1) 

o2 

PRVR(k  -1) 

+ 

ap/T 

[QrrW] 
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%(*) 

(1  -  a)2  2(1  -  a)2T  (1  -  a)2T2 

pdve(k) 

- 

-0(1  -  a)/T  (1  -  a)(l  -  20)  (l-a)(l-0)T 

pvov0(k) 

( 0IT )2  -20(1- 0)1  r  ( 1-0 )2 

Pgeik- 1) 

a2 

X 

peve(k-l) 

+ 

U0/T 

pveve(k~D 

( 0IT )2 
—  _ 

[Qs0(fe)] 


aRn  =  PRR<k)  +  2U!m)TPRVR(k)  +  (j/m )2T2  P vrVr (k) 


<j|d  =  Pee(k)  *  2{jlm)TP6Vg(k)  +  Ulm)*T* PVgVg{k) 


2rp2 


PRpQp  =  0  because  8m  and  Rm  are  uncorrelated. 


(4.23) 

(4.24) 

(4.25) 

(4.26) 


Since  Eqs.  (4.3),  (4.4),  (4.6),  and  (4.7)  are  identical  with  Eqs.  (4.22)-(4.25)  in  form. 


°R  ~ 
Rp 

Fa2 

(4.27) 

li 

nq?1 

D 

Fop 

°m 

(4.28) 

PRfflp  ~ 

0 

(4.29) 

Pm 

(4.30) 

= 

em . 

(4.31) 

Eqs.  (4.13H4.17)  are  identical  to  Eqs.  (4.27)-(4.3l),  a  fact  which  shows  that  the  system 
of  Fig.  2.1  is  the  same  as  that  of  Fig.  4.1.  To  justify  this  result  the  following  argument 
is  given.  The  system  shown  in  Fig.  2.1  is  redrawn  in  Fig.  4.2.  Recall  that  in  Section  2.0 
the  polar-to-cartesian  and  cartesian-to-polar  transformations  were  shown  to  be  approxi¬ 
mately  linear  over  the  region  governing  the  means  and  covariances.  The  order  of  operation 


Fig.  4.2— Filter  system 
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for  linear  operators  can  be  changed,  yielding  the  system  shown  in  Fig.  4.3.  For  stationary 
target  case  the  two  transforms  caned,  yidding  Fig.  4.1. 


Fig.  4.3— Interchange  of  linear  operations 


The  value  of  F  can  be  computed  as  follows:  P(A)  is  set  equal  to  P(fe  - 1)  and  the 
resulting  algebraic  equations  are  solved; 


PR^k) 

pee(k) 

2 P  -  3ap  +  2a2 

(4.32) 

4 

On 

"m 

a(4  -  2a-  p) 

PRVR(k) 

_  P6Ve(V 

P(2  a-p) 

(4.33) 

°R 

aB 

“m 

a(4-  2a -p) 

pvRvR(k)  pVsVg(k )  012a2 -a3  +  2/J-o$] 


UR  r, 


On, 


oc(4-2a-P) 


(4.34) 


F 


PrrW 


+  (2//m) 


PRVR(k) 


U.'my 


PVRVR(k ) 


Rr 


(4.35) 


Hie  peak  variance  occurs  for  j  —  m.  In  Fig.  4.4,  F  is  plotted  vs  a  and  ft  for  j  =  m.  Ob¬ 
serving  this  figure  one  finds  that  the  amount  of  smoothing  (a  and  0 )  controls  the  peak 
noise  levds.  In  addition,  the  noise  varies  between,  each  scan  of  the  search  radar  from  k 
to  k  =  1.  This  is  plotted  vs  j  in  Fig.  4.5  for  large  m.  Observing  this  figure  one  finds  that 
the  noise  is  nonstationary  but  periodic  with  time  under  steady-state  conditions. 


4.3  Some  Nonstationary  Solutions 

A  target  is  flown  at  a  constant  vdocity  in  a  straight-line  trajectory  as  shown  in  Fig. 
4.6.  The  measurement  standard  deviations  are  assumed  to  be  oRfn  =  250  ft  and  o$m  = 
0.5°.  Equations  (2.10H2.14)  are  used  to  obtain  o|m,  Pxmymi  in  terms  of 

the  measurement  variances  and  target  trajectory.  Using  these  results  with  Eqs.  (4.1), 
(4.2),  and  (4.3H4.8),  we  find  aRp,  o£p,  pRjfip,  P.p,  and  Qp.  The  results  for  several  trajec¬ 
tories  and  values  of  smoothing  coefficients  are  shown  in  Figs.  4.7-4.12.  In  all  cases  only 
the  envelope  of  the  peak  noise  (/=  m)  is  shown.  The  covariances  will  have  a  ripple 
between  samples  (k)  and  (k  + 1).  In  addition  the  dotted  lines  show  the  covariances  if  the 
target  is  stationary.  Figures  4.7-4.12  show  that  the  output  noise  processes  are  in  general 
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Fig.  4.4— Predicted  position  noise  power  as  a 
function  of  (a, 0) 


Fig.  4.5— Normalized  predicted  azimuth  variance 
as  a  function  of  time 
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Fig.  4.8— Covariance  of  predicted  range 
ard  azimjuth  for  straight-line  target  tra¬ 
jectory  Xm  =  2  n.mi.,  v  -  2000  ft/s, 
a  =  0.1,  0  =  0.005 


J _ I _ ill _ ! _ L 


-90  -45  0  45  90 

Bp  (Degrees) 

Fig.  4.9— Covariance  of  predicted  range 
and  azimuth  f  jt  straight-line  trajectory 
Xm  =  2  r.mi.,  v  =  2000  ft/s,  T  =  8  s, 
a  =  0.56,  $  =  0.85 


nonstationary  and  depend  upon  the  target  trajectory  and  velocity,  the  sampling  time,  and 
the  filter  parameters. 

A  circular  flight  path  is  flown  as  shown  in  Fig.  3.2.  The  covariances  of  the  predicted 
range  and  azimuth  are  shown  in  Figs.  4.13-4.16  as  a  function  of  time  for  various  condi¬ 
tions.  Again  only  the  envelope  of  the  peak  variances  is  plotted  Kj=  m).  At  the  far 
ranges  or  low  velocities  the  variances  of  range  and  azimuth  approach  the  stationary  solu¬ 
tion  values,  although  at  the  long  ranges  it  was  found  that  the  correlation  did  not  go  to 
zero  but  was  a  function  of  the  turning  motion.  Again  it  is  found  that  the  covariances 
can  be  a  function  of  target  trajectory,  velocity,  and  range;  sampling  time;  measurement 
uncertainty;  and  filter  parameters. 

The  effect  of  the  faster  sampler  {/  =  1, . . .  m)  is  shown  in  Fig.  4.17.  Observing  this 
figure  one  finds  that  the  covariances  ripple  between  the  scan  time  of  the  search  radar  in 
a  similar  manner  as  found  for  the  stationary  solutions. 
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Fig.  4.10— Covariance  of  predicted  tangs 
and  azimuth  for  straight-line  trajectory 
Xm  =  30  n.mi.,  v  =  2000  ft/s,  a  =  0.56, 
P  =  0.85 


Fig.  4.11— Covariance  of  predicted  range 
and  azimuth  for  straight-line  target  tra¬ 
jectory  XK  =  30  n.mi.,  v  =  500  ft/s, 
a  =  1.0,  P  =  1.0 


4.4  Discussion  of  Results 

The  covariance  equations  were  found  for  both  the  polar  and  cartesian  coordinate 
filters.  A  simple  closed-form  solution  was  found  for  the  covariances  in  both  filters  when 
the  target  was  stationary  in  space.  For  these  cases  the  output  correlation  was  zero  and 
the  output  variances  only  depended  upon  a,  P,  and  the  input  measurement  variances. 

For  nonstationary  targets  the  cartesian  coordinate  filter  yielded  output  covariances  which 
depended  upon  a,  P,  input  measurement  variances,  target  trajectory  and  speed,  and 
sampling  time.  In  general  the  covariances  increased  as  one  came  near  the  radar  and 
tended  to  approach  the  statioi,  ary  target  solutions  at  the  far  ranges  except  that  the  cor¬ 
relation  depended  upon  the  turning  motion. 


5.0  COMPARISON  OF  POLAR  AND  CARTESIAN  COORDINATE  a-0  FILTERS 
OPERATING  UNDER  SHORT  FADE  CONDITIONS 

This  section  is  concerned  with  evaluating  the  tracking  performance  of  the  polar  and 
cartesian  a-P  filters  for  short  fade  conditions  and  for  the  track  handoff  problem.  Section 
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Fig.  4,12— Covariance  of  predicted  range  and  azimuth  for 
straight-line  trajectory  flying  directly  at  radar:  v  *  2000 
ft/s,  tt  =  0.56,  /?  *  0.85 
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Fig.  4.13— Covariance  of  predicted  range  and 
azimuth Qfor  circular  flight  path  Rm  =  21.2  n.mi„ 
$n  =  45°,  »  =  2000  ft/s,  3-g  turn,  a  =  0.56,0  =  0.85 


lr- 


1 


'0 


j 


80 


160 

TIME  (S) 


240 


320 


NRL  REPORT  7S48 


31 


t 


S 


Fig.  4.16— Covariances  of  predicted  range  and 
azimuth  for  circular  flight  path  Rm  =  21.2  n.mi., 
-  45°,  v  =  2000  ft  Is,  3-g  turn,  a  =  0.1 ,0  =  0.005 
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Fig.  4.17— Covariances  of  predicted  ranjt  and 
azimuth  for  circular  flight  path  including  inter¬ 
sample  ripple:  Rm  =  21.2  njnl,  =  45°,  v  = 
2000  ft/s,  3-g  turn,  a  «*  0.56,  0  =  0.85 
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5.1  describes  the  mean  errors  due  to  short  fades,  and  Section  5.2  describes  the  effects  on 
the  covariances.  Section  5.3  presents  the  results  of  the  simulation  on  the  track  handoff 
problem. 

5.1  Mean  Errors 

The  predicted  value  strategy  used  for  processing  information  under  fading  conditions 
is  defined  as  follows.  When  a  fade  occurs  the  predicted  position  is  set  equal  to  the  meas¬ 
ured  position.  The  a-@  filter  equations  in  cartesian  coordinates  given  in  Eqs.  (3.1H3.4) 
reduce  to 


~X(k)~ 
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T 

X(k-1) 

yx(k) 
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yx(k- 1) 

(5.1) 


Xp(k  +  j/m)  =  [1  .  jTIm) 
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and 


Yp(k  +  ilm)=ll  jTfM) 


Y(k) 

Vy(k) 
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during  the  period  of  the  fade.  Otherwise  they  are  the  same  as  before.  Observing  Eqs. 
(5.1)-(5.4),  one  finds  that  they  are  of  the  same  form  as  Eqs.  (3.1H3.4)  except  that 
a  =  0  =  0.  Therefore  Eqs.  (3.1)-(3.4)  can  be  used  to  represent  the  system  under  fading 
conditions,  except  that  a  and  0  are  time  varying  between  these  set  values  and  a  =  j3  =  0. 

In  a  similar  manner  the  filter  in  polar  coordinates  is  found  to  be  represented  by  Eqs. 
(3.16H3.19)  with  time-varying  coefficients  of  a  =  0  =  0  and  the  original  set  values  of  a 
and  0  corresponding  to  a  fade  and  no-fade  condition.  Under  fading  conditions  the  polar 
coordinate  filter  moves  the  target  along  curvilinear  lines  at  the  last  known  radial  and 
traverse  velocities  during  the  fade.  The  cartesian  coordinate  filter  moves  the  target  along 
straight  lines  at  the  last  known  velocity. 

Several  constant-velocity,  straight-iine  target  trajectories  as  shown  in  Fig.  4.6  were 
flown.  In  all  cases  using  the  cartesian  coordinate  filter,  the  mean  error  between  the  pre¬ 
dicted  and  true  target  positions  was  zero  under  the  fading  sequences.  By  using  the  polar 
coordinate  filter,  the  mean  error  between  the  predicted  and  true  target  azimuths  was 
found  under  several  conditions  as  shown  in  Figs.  (5.1H5.4).  A  one  in  the  fading  sequence 
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Fig.  5.1— Error  between  predicted  and  true  target  positions  at  sampling 
instants  k  (j  =  m)  for  straight-line  target  trajectory  under  fading  conditions 
(tracking  in  polar  coordinates} 
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Fig.  5.2— Error  between  predicted  and  trun  ••srgel  positions  at  sampling 
instants  k  (j  =  m)  for  straight-line  target  trajectory  under  fading  conditions 
(tracking  in  polar  coordinates) 
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Fig.  5.3— Error  between  predicted  and  true  target  positions  at  sampling 
instants  k  (j  =  m)  for  straight-line  target  trajectory  under  fading  conditions 
(tracking  in  polar  coordinates} 


FI’  5.4— Errors  between  predicted  and  true  target  position  at  sampling 
instant"  k  (j  =  m)  for  straight-line  target  trajectory  under  fading  conditions 
(tracking  in  polar  coordinates) 
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represents  a  fade,  whereas  a  zero  corresponds  to  no  fade.  In  all  cases  the  envelope  of  the 
errors  at  sample  instants  (;'  =  m)  is  shown,  ignoring  the  ripple  errors  due  to  the  faster 
sampler.  These  figures  indicate  that  the  error  becomes  larger  as  the  target  comes  closer 
to  the  radar  and  when  the  filter  uses  heavier  smoothing.  In  addition  the  error  grows  dur¬ 
ing  the  time  the  fade  is  present.  The  reason  this  occurs  is  that  the  target  is  being  pro¬ 
jected  along  curvilinear  lines  during  the  fades  and  is  moving  in  a  straight  line.  The  error 
situation  would  be  reversed  if  the  target  were  moving  along  a  traverse  line  rather  than  a 
straight  line. 

A  target  is  flown  in  a  circular  trajectory  and  the  mean  errors  between  the  predicted 
and  true  positions  are  computed  for  a  given  situation  including  a  sequence  of  fades.  The 
results  are  shown  in  Figs.  5.5  and  5.6.  These  figures  show  that  at  least  at  the  farther 
ranges  and  short  fading  conditions  the  error  is  approximately  the  same  using  either  the 
polar  or  cartesian  coordinate  system  filter. 

In  general  it  appears  that  at  the  farther  ranges  where  the  near  effect  accelerations 
are  not  present  in  the  polar  coordinate  system,  the  mean  errors  in  either  filter  system 
under  short  fading  conditions  are  nearly  the  same.  The  polar  coordinate  filter  performs 
better  for  constant-velocity  targets  moving  in  the  polar  directions,  whereas  the  cartesian 
coordinate  filter  performs  better  for  targets  moving  in  straight  lines. 


Fig.  5.5— Mean  error  between  predicted  and  true  positions  for  a  circular 
target  trajectory  using  an  a-0  filter  in  cartesian  coordinates  under  fading 
"renditions 
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Fig.  5.6— Mean  error  between  predicted  and  true  positions  for  a  circular 
target  trajectory  using  an  a-0  filter  in  polar  coordinates  under  fading 
conditions 


5.2  Covariance  Description 

Under  fading  conditions  and  using  the  predicted  position  strategy  as  outlined  in  Sec¬ 
tion  5.1,  tha  covariances  are  described  in  the  same  manner  as  in  Section  4.0  except  that 
when  a  fade  occurs  a  and  are  set  equal  to  zero.  A  constant-velocity  target  is  flown  in  a 
straight  line  and  the  covariances  are  computed  under  a  fading  condition  for  two  cases  as 
shown  in  Figs.  5.7  and  5.8.  The  cartesian  coordinate  filter  is  used  and  the  envelope  of 
the  covariances  (j  =  m)  is  shown.  The  covariances  increase  during  a  fading  condition. 

The  eonvariances  were  computed  using  the  polar  coordinate  filter  and  are  shown  in  Fig. 
5.9.  Unlike  the  cartesian  coordinate  filter,  in  this  case  the  covariances  are  independent  of 
target  trajectory,  speed,  and  sampling  time. 

A  circular  flight  path  is  flown  and  the  covariances  are  computed  using  the  polar  and 
cartesian  coordinate  filters.  The  results  are  shown  in  Figs.  5.10-5.12,  which  show  that 
the  covariances  increase  as  the  target  fades. 
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Fig.  5.7— Covariances  of  predicted  range  and  azimuth  for 
straight-line  target  trajectory  Xm  -  30  njni.,  u  -  2000  ft/s, 
T-  4,  a  =  0.56,0  =  0.85,  for  fade  sequence  (10001000100....) 
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Fig.  5.8— Covariances  of  predicted  range  and  azimuth  for 
straight-line  target  trajectory  Xm  =  2  njna,  v  =  2000  ft/s,  T  = 
4,  a  =  0.56,  0  =  0.85,  for  fade  sequence  (10001000100 _ ) 
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Fig.  5.9— Covariances  of  predicted  range  and  azimuth  for  a 
filter  in  polar  coordinates  operating  under  fading  conditions 
(1000100010...),  a  =  0.56,  j?  =  0.85 


Fig.  5.10— Covariances  of  predicted  range  and 
azimuth  for  circular  flight  path  Rm  =  21.2 
n.mi.,  0m  =  45°,  v  =  2000  ft/s,  T  =  4  s,  3-g 
turn,  a  =  0.56,  0  =  0.85,  under  fading  condi¬ 
tions  (10001000100...),  (X-Y)  track 
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Fig.  5.11— Covariances  of  predicted  range  and 
azimuth  for  circular  flight  path  Rm  =  210.2 
n.mL,  =  45°,  v  =  2000  ft/s,  T  =  4  s,  3-g 
turn,  a  x  0.56,  P  =  0.85,  under  fading  condi¬ 
tions  (1000100010...),  (X-y)  track 


12 

I-  *t>AAAAAAA AAJ 

0I - 1 - JL. - 1 - 1 

0  80  160  240  320 

TIME  (Si 


12 

>■ '  -AAAAAAAAA/ 

0i - 1 - l - 1 - 1 

0  80  160  240  320 
TIME  (S) 


0  60  160  240  320 

TIME  (S) 


TIME  (S) 


TIME  (S) 


Fig.  5.12— Covariances  of  predicted  range  and 
azimuth  for  a  filter  in  polar  coordinates  oper¬ 
ating  under  fading  conditions  (10010010...) , 
a  =  0.56,  p  =  0.85 
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5.3  Track  Handoff  Problem 

In  this  section,  the  abilities  of  the  polar  and  cartesian  a-P  filters  to  perform  the  track 
handoff  problem  are  compared.  Since  the  details  of  the  track  handoff  simulation  are 
given  in  Ref.  2,  only  the  basic  facts  will  be  presented  here. 

The  geometry  of  the  situation  is  shown  in  Fig.  5.13.  The  coordinates  of  the  target 
at  time  f  =  0  are  (xq,  y0,  fig),  its  height  when  crossing  the  y  axis  is  hf,  its  ground  speed  is 
v,  and  its  heading  is  A.  The  radar  coordinates  are  (0, 0,  hr). 


TARGET 


Fig.  5.13— Geometry  of  radar  and  target 


The  simulation  is  run  in  the  following  manner:  The  target  is  assumed  to  be  detected 
on  every  scan  of  the  search  radar.  Initially,  three  target  positions  are  generated,  and  the 
a-p  filter  is  initialized.  After  each  additional  detection,  the  a-P  filter  is  updated  and  is 
used  to  continuously  estimate  the  target’s  coordinates.  Starting  with  the  third  sample,  the 
center  of  the  tracking  scan  pattern  is  centered  on  the  predicted  position  of  the  target. 

The  search  pattern  of  the  tracking  radar  is  initialized,  and  during  each  update  time,  the 
program  calculates  whether  or  not  the  target  is  located  within  the  acceptance  beam  of  the 
tracking  radar  (Section  5.4).  The  simulation  continues  until  the  target  crosses  the  y  axis. 
The  output  of  a  single  case  is  a  series  of  correct  and  incorrect  handoffs  between  the 
search  and  tracking  radars.  Many  cases  are  runs,  and  the  probability  of  handing  off  as  a 
function  of  range  is  estimated. 

The  initial  simulation  was  run  with  the  following  target  parameters:  x0  was  uni¬ 
formly  distributed  between  121,600  ft  (20  n.mi.)  and  122,600  ft,  y0  =  0,  h0  =  10,000 
ft,  hf  =  5000  ft,  A  is  uniformly  distributed  between  5.9°  and  6.1°,*  and  V  is  uniformly 
distributed  between  2100  and  2300  ft/s.  The  radar  is  at  a  height  of  80  ft,  has  a  range 
resolution  of  250  ft  and  a  scanning  rate  (update  time)  of  4  s,  and  measures  the  azimuth 
position  with  a  standard  deviation  of  0.5°.  The  deviation  accuracy  of  the  radar  will  be 
varied  in  the  simulation.  It  will  have  a  standard  deviation  of  1°,  or  else  the  radar  (a  2-D 
radar)  assumes  that  the  elevation  of  the  target  is  always  12°  (the  bottom  of  the  24°  scan 
pattern  of  the  tracking  radar  is  set  on  the  horizon).  The  filter  parameters  are  a  =  0.6 
and  P  =  0-9  for  range,  a  =  ft  =  0.5  for  azimuth,  and  a  =  0.56  and  p  -  0.85  for  X  and  Y. 


*This  makes  the  target  pass  within  2  mi  of  the  radar. 
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For  each  of  the  two  elevation  accuracies  and  two  filters,  50  cases  were  run;  and  the 
probability  of  the  target  being  in  the  beam  of  the  tracking  radar  on  the  last  scan  pattern, 
vs  the  target  range,  is  shown  in  Figs.  5.14  and  5.15.  It  is  obvious  that  one  can  hand  off 
targets  at  closer  ranges  using  the  X-Y  filter.  This  is  because  oi  the  large  accelerations  in 
the  R-Q  coordinate  system  for  crossing  targets. 
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Fig.  5.15— Probability  of  handoff  using  X-Y  filter 
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5.4  Discussion  of  Results 

The  mean  errors  in  the  polar  and  cartesian  filters  were  found  to  increase  during  fades 
unless  the  target  was  moving  at  a  constant  velocity  along  the  polar  or  cartesian  coordi¬ 
nates,  respectively.  For  short  fades  at  the  ionger  ranges,  there  seemed  to  be  very  little 
difference  between  the  mean  errors  in  either  system.  The  covariances  increased  during 
fading  conditions. 

One  can  hand  off  targets  at  closer  ranges  using  the  X-Y  filter.  This  is  because  of  the 
large  accelerations  in  the  R-d  coordinate  system  for  crossing  targets. 


6.0  CONCLUSION 

This  report  describes  an  approximate  analytical  procedure  for  determining  the  errors  in 
an  a-p  filter  operating  in  cartesian  coordinates.  These  errors  are  separated  into  mean  and 
covariance  errors  and  a r*  compared  to  the  errors  in  an  a-p  filter  operating  in  polar  coordinates. 

The  polar  to  cartesian  and  cartesian  to  polar  coordinate  transformations  are  shown 
to  be  approximately  linear  over  the  space  in  which  the  central  and  near  tail  regions  of 
the  probability  density  lie.  Since  the  measurement  probability  density  is  Gaussian,  it  is 
shown  that  the  probability  densities  after  the  transformations  can  be  approximated  with 
a  high  degree  of  accuracy  with  Gaussian  densities  as  long  as  the  far  tail  region  is  of  little 
concern.  Using  this  result,  we  find  the  covariances  at  the  output  of  the  filters.  Closed- 
form  steady-state  solutions  are  found  for  the  covariances  in  the  polar  coordinate  filter 
and  for  stationary  targets  using  the  cartesian  coordinate  filter.  These  covariances  depend 
upon  a,  P,  and  measurement  variances.  For  moving  targets,  the  cartesian  coordinate  filter 
yields  output  covariances  which  are  nonstationary.  Their  values  depend  upon  a,  P,  meas¬ 
urement  variances,  target  trajectory,  target  speed,  and  sampling  time.  In  steady  state  and 
for  the  far  ranges,  the  output  variances  using  the  cartesian  coordinate  filter  approach  the 
variances  obtained  from  the  polar  coordinate  filter.  At  the  near  ranges,  the  covariances 
are  in  general  larger  in  the  cartesian  coordinate  filter  as  compared  to  the  polar  coordinate 
filter. 

At  the  far  ranges  the  mean  errors  using  either  filtering  system  are  essentially  the 
same.  But  at  the  close  ranges  the  mean  error  in  the  polar  coordinate  system  is  in  general 
larger  than  in  the  cartesian  coordinate  filter. 

In  the  study  on  the  effects  of  short  fades  using  the  predicted  target  strategy,  it  was 
generally  found  that  both  the  mean  errors  and  covariances  increased  during  the  time  of 
the  fade. 

In  the  study  of  the  track  handoff  problem  it  was  found  that  one  could  hand  off 
targets  at  closer  ranges  using  the  X-Y  filter.  This  is  because  of  the  large  accelerations  in 
the  R-P  coordinate  system  for  crossing  targets. 
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Appendix 

CALCULATION  OF  MEANS  AND  COVARIANCES  OF  RADAR 
MEASUREMENTS  IN  TERMS  OF  CARTESIAN  COORDINATES 


A  radar  measurement  is  given  in  range  R  and  azimuth  0,  where  R  and  0  are  assumed 
to  be  uncorrelated,  Gaussian-distributed  random  variables  with  means  0  and  R,  and 
variances  o£  and  o§.  The  problem  is  to  determine  the  means  and  covariances  of  the 
two  quantities  X  and  Y  defined  as 


R  cos  6 

<A1) 

R  sin  0  . 

(A2) 

In  the  calculations  the  following  two  facts  are  used  extensively.  The  approximation 


(-2o|)N 

Nl 


(A3) 


is  used  because  the  azimuth  standard  deviation  Og  in  radians  found  in  typical  search  radars 
is  small.  The  integral 

a 2 

— ,  where  a  >  0  (A4) 


r 


g-a2*2  CQS  bx  dX  = 


V?  e"*2'4 
2a 


is  used  in  each  of  the  calculations. 

The  major  steps  in.  the  calculations  are  the  following. 

Mean  of  X,  X 

X  =  E[R]  E[ cos  0]  where  E[R ]  =  R 

E{ cosSj  - - - -  f  cos 9  <f  (1/2)[(0-^)/as)2]  dd  . 

y/M  Oq  -Loo 

Change  of  variable  «  =  9  -  9 

IS  [cos  0]  - - - —  f  [cos  to  cos  6  —  sin  co  sin  9]  e(“U2)(to/Ofl) 

y/Z*  Oq  J_oo 
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E [cos  01  =  2  C- —  I  cos  we 


r 


(i/VTo^w2  dw 


V2?  og  jo 

E[ccs0]  =  e_a«/2  cos  0 

X  -  R  e-a®  ^2  cos  0  ;  for  small  Og  ,  X  =  R  cos  0 
Mean  of  Y,  Y 

Similar  to  calculation  of  X 

Y  =  R  e~°®  sin  0  ;  for  small  og  ,  Y  =  E  sin  0 
Variance  of  X,  VAR  X 

VarX  =  £[X21  -  (£[X1)2  =  E[R2)  E[ cos2  0]  -  X2 
iJ[cos*  «,  =  f  cos*  9  e-(l/2)l(»-f)/o.)2l  d9 

yJW  Og  loo 

£tC0S2  „  i  x  _i_  r  m 

y/ 2tT  Oq  'Loo 


ax  i_  r 

2  V2tF  og  loo 


cos 


20  e-(l/2)[(0-5)M>)2]  d0  . 


Let  w  =  0  -  0  ; 


1  ,  1 


£ [cos^  01  =  -5-+ 


J  cos(2w  +  20)e(  1'2Hw/0®)  dco 


2  2  V2Fog  1 

„00 

£  [cos2  01  =  —  +  -i-  X  - - -  I  (cos  2w  cos  20  -  sin  2w  sin  20) 

2  *  y/ZS  Og  loo 

x  e(-l/2)(w/o„)2  d(ji 


Elcos^  01  =  2  + 


1  2  cos  20 


—  ao 


2  V2ir  og  Ico 


L 


2w  e-U/^V2)2«2  d£J 


cos 


Eq.  (2.13) 


Eq.  (2.14) 


46 


BEN  a  CANTBELL 


£[cos2  0]  =  i  +  \  e~2a°  cos  20 
E[R2]  =  4  +  (R)2 

Var  X  =  £[Jt2]  £[cos2  0]  -  (X)2  ;  for  small  %  , 

Var  X  =  o|  cos2  0  +  (R)2  c|  sin2  0 
Variance  of  Y,  Var  Y 

Similar  to  calculation  of  Var  X,  for  small  Og  , 

Var  Y  =  o|  sin2  0  +  (R)Z  of  cos20  . 

Covariance  of  X  and  Y,  Coo  XY 

covXY  =  E[XY]  -  E[X]  jB[Y]  =  E[R2]  £[cos0  sin0]  -  XY 
E[ cos 0  sin  0)  *  4  x  — - -  f  sin  20  e-(1/2)[(0-^)/CT9)2]  d6 

z  0*)  loo 

CO  =  0  -  V 

<30 

£[cos0  sin0]  =  — —  f  (sin  2a>  cos  20 
2  %/2?  ag  -ioo 

f  cos  2 to  sin  20)  dco 

£[cos0  sin0]  =  — f  cos  2co  e  ^  ‘  '  Y  ^  w2  dco 
'/2w  o9  Joo 

J5[cos0  sin  0  ]  =  -|-  e~2ofl  sin  20 
EIR2]  =  a2  +  (fi)2 

covXY  =  £(J?21  £[cos0  sin0]  -  XY;  for  small  og  , 
o2  So2 

cov  XY  -  —•  sin  20 - sin  20  . 


Eq.  (2.10) 


Eq.  (2.11) 


Eq.  (2.12) 


